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I.  INTRODUCTION 


The  aerodynamic  performance  of  flared  projectiles  has  not  been 
investigated  in  detail  -  either  by  the  wind  tunnel  techniques,  which  usually 
focus  on  the  effects  of  perturbations  of  the  flare  geometry;  or  on  the  free 
flight  ranges,  which  normally  employ  the  flare  for  its  strong  stabilizing 
influence  in  order  to  survey  some  other  performance  parameter.  For  a  given 
projectile,  the  wind  tunnel  experiments  can  efficiently  determine  axial  force, 
normal  force  and  spin  deceleration  characteristics  and  the  influence  of  the 
Reynolds  number  on  these  coefficients.  However,  the  dynamic  coefficients, 
and  CM„  +  C  __  ,  are  not  easily  established  and  the  axial  force 

coefficient  is  clearly  not  equivalent  to  the  drag  coefficient  which  is 
required  for  application  design.  Free  flight  ranges  are  most  suitable  for 
measurements  of  the  drag  and  the  dynamic  coefficients  and,  ultimately,  the 
dispersion  and  accuracy. 

The  few  flare  control  data  available  verify  that  the  high  drag  associated 
with  the  flare  makes  such  a  projectile  unattractive  for  anything  but  very 
short  range  applications,  although  a  recent  publication  in  the  form  of  a 
careful  computational  study  describes  a  very  novel  flare  concept  which  will 
require  evaluation  by  experimental  methods. 

This  present  report  is  intended  to  be  used  in  the  formulation  of  quick 
estimates  of  the  static  aerodynamic  characteristics  of  the  conventional  flared 
projectile.  The  data  is  based  primarily  on  References  2  and  5  and  the 
exposition  follows  that  of  the  corresponding  finned  projectile 
presentation.  The  drag,  static  moment  and  normal  force  coefficients  are 
determined  from  the  simplified  geometry  of  the  projectile  and  the  bounding 
aerodynamic  conditions,  i.e.,  sea  level,  flat  fire,  2<M<5,  and  small  yaw.  The 
computer  program  is  written  in  the  Hewlett-Packard  Basic  language  and  is 
easily  convertible  to  other  Basic-language  microcomputers. 

II.  PROCEDURE  AND  RESULTS 

The  empirically-developed  expressions  for  the  aerodynamic  contributions 
of  the  components  are  based  on  those  of  the  corresponding  elements  for  the 
similar  finned  projectile  -  following  the  form  of  Reference  5.  Modifications 
are  introduced  as  indicated  by  the  test  data  of  Reference  2.  Caliber 
notation,  where  a  representative  linear  dimension  indicates  a  reference  length 
and  the  mass/force  dimension  is  converted  to  a  specific  gravity  equivalent,  is 
used  in  the  physical  description  and  in  the  interior  of  the  program  operation. 

Figure  1  is  a  sketch  and  Figure  2  is  a  set  of  shadowgraphs  of  a  flared 
projectile.  Tables  1,  2  and  3  define  the  equations  used  to  describe  the  aero¬ 
dynamic  parameters.  The  initialization  instructions  are  presented  in  Appendix 
A  and  the  full  program  listings  follow.  Figures  3. a  through  3.g  present  a 
projectile  outline,  input  tabulation,  a  table  of  the  static  aerodynamic  co¬ 
efficients  and  graphs  of  the  normal  force  slope,  drag,  static  moment  slope, 
lift  force  slope,  accuracy  factor  and  velocity  decrement.  The  range  is  2<M<5. 
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TABLE  2.  NORMAL  FORCE  AND  STATIC  MOMENT  SLOPE  COEFFICIENT 


Normal  Force  and  Static  Moment  Coefficients 


ft  —  7 

Normal  Force  C^g  -  (l.9  +  1.3-j-  +  .0149  ~)[6  >'.6103 

Coefficient  ‘  n 

O  1  ?  e  |  Q  Q 

Center  of  Pressure  x  .  -  (.89  +  .75  —  +  .5  — )  M""  6" ‘ 

c.p.B  1  S  ; 

of  Normal  Force  n 


Static  Moment  C^  -  (*c  „)(C  B) 


Coefficient 


Normal  Force 
Coefficient 


cNaF  '  ^2>67  "  d  ■ •  C«o«e)CM/B) 


Center  of  Pressure  x  „  -  £  +  £  +  .1  (d.  -  I.’/tanB 

c.p.F  n  a  1  b  ' 

of  Normal  Force 


Static  Moment 
Coefficient 


Normal  Force 
Coefficient 

Static  Moment 
Coefficient 


W  *  ^Xc.p.F^CNoF^ 


So  ’  CNoB  +  CNoF 


Sfa  ”  SfciB  +  SfaF 


Center  of  Pressure  x  -  C_  / C.. 

c .  l  .  flo  Nc 

of  Normal  Force 


Normal  Force 
Coefficient 


C  *  C 
No  No 


Center  of  Pressure  x 

c 

of  Normal  Force 


c.p.  c.p. 


Static  Margin  I  ix  ■  x 


c.p.  c.g. 


Static  Moment 


C.  -  (x  -  x  )  (C  ) 
no  c.p.  c.g.  No 


TABLE  3.  RETARDATION  AND  ACCURACY  FACTOR 


Retardation 


Mach  Number 
along  Trajectory 

M  -  b 

1  R  e  Q8  -  c 

Average  Velocity  Decrement 

iv  -  M°  "  ^  (V)  . 

s  sonic 

Accuracy  Factor 


5  CMO  «i2 


Accuracy  Factor 


ft 


.967' 

2.068 


4 


3 . 247p2Q0^5  19 


15  DEG  FLRRE 


NT  =  4.11  WCRL3 
IX  -  0  ICRL5 
IY  -  2.08  ICRL5 
DIfl-  2.757  IN/CRL 


ALL  VALUES  ARE  IN  CALIBERS  UNLESS  OTHERWISE  NOTED 

CONICAL  NOSE  LENGTH:  .96? 

CYLINDRICAL  BODY  LENGTH:  1.101 

GROOVE  LOCATION /NOSE:  0 
GROOVE  LENGTH:  0 
MAX  FLARE  DIAMETER:  1.632 
FLARE  ANGLE:  15 
0.00 

CENTER  OF  GRAVITY:  1.758 

MACH  NUMBER  AT  MUZZLE:  3 

MAXIMUM  RANGE  (METERS  > !  1000 

PROJECTILE  WEIGHT  <CAL3>:  4.11 

AXIAL  MOMENT  OF  INERTIA  <CAL5>:  0 

TRANSVERSE  MOMENT  OF  INERTIA  (CALS):  2.08 

1.0  CALIBER  REFERENCE  DIAMETER  (IN.  2.757 


Figure  3a.  Computer-Generated  Output:  Projectile  Outline  and  Input  Parameters 


CNAB 

2.0 

STATIC  AERODYNAMIC  COEFFICIENTS  FOR 
LONG  ROD  FLARED  PROJECTILES 

MACH  NUMBER 

2.5  3.0  3.5  4.0 

4.5 

5.0 

1.976 

1.911 

1.888 

1 . 383 

1.887 

1 . 898 

1.911 

XCPB 

1.331 

1 . 507 

1.631 

1 . 730 

1.817 

1 . 394 

1 . 966 

CMAB 

2.629 

2.880 

3.079 

3.253 

3.428 

3.594 

3.757 

CNAF 

2.520 

2.381 

2.315 

2.277 

2.254 

2.238 

2.227 

CNRT 

4.496 

4.293 

4.203 

4.  160 

4.  141 

4.  136 

4.139 

CMRT 

3.435 

8. 366 

3.412 

8.505 

8.621 

8.751 

8.889 

CG-CP 

-.  118 

-.  191 

-.244 

-.286 

-.324 

-.358 

-.390 

CMRT , CG 

-.532 

-.819 

-1.023 

-1.191 

-1.341 

-1.480 

-1.613 

CDWN 

.611 

.  574 

.545 

.  522 

.503 

.487 

.  473 

CDB 

.527 

.452 

.377 

.  302 

.227 

.  153 

.  078 

CDV 

.046 

.  041 

.037 

.032 

.  027 

.  022 

.013 

CDGR 

0.000 

0 . 000 

0.000 

0.000 

0.000 

0.000 

0.000 

CDWF 

.  150 

.  120 

.  100 

.  086 

.  075 

.  067 

.  860 

CDT 

1 . 334 

1 .  187 

1.059 

.  942 

.  833 

.  729 

.  628 

CLA 

3.  162 

3.105 

3.  143 

3.218 

3.303 

3.407 

3.510 

J  ZETfl 

-3.009 

-1.918 

-1.554 

-1.367 

-1.249 

-1 . 165 

-1. 102 

•®i 
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Figure  3b.  Computer-Generated  Output:  Table  of  Aerodynamic  Coefficients 


MRCH  NUMBER 


Figure  3c,  Computer-Generated  Output:  Normal  Force  Slope  Coefficient 
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MRCH  NUMBER 


Figure  3e.  Coaputer-Generated  Output:  Static  Moaent  Slope  Coeffici< 
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MUZZLE  MRCH  No.  3.000,  INITIAL  VAN  CYCLE  1 1 . 509M 


2000 


1750 


1500 


if)  1250 
\ 

2: 

1 000 

u  750 
> 


RANGE  (KM) 

RANGE  0  M ,  RETARD AT I ON  13  40.284  M/S/KM 
RANGE  500  M,  RETARDATION  1 340 . 284  M/S/KM 
RANGE  1000  MP  RETARDATION  927.116  M/S/KM 
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APPENDIX  A 


INITIALIZATION  INSTRUCTIONS 

To  use  this  program,  the  operator  requires  only  the  entry  of  the  geometry 
and  the  physical  properties  of  the  projectile,  the  range  excursion  and  the 
muzzle  Mach  number.  Caliber  notation,  where  a  representative  linear  dimension 
provides  a  reference  length  and  the  mass /force  dimension  is  converted  to  a 
specific  gravity  equivalent,  is  employed  extensively.  The  reference  diameter 
may  be  in  inches  or  in  millimeters  and  the  range  input  in  meters. 

The  program  will  call  for  the  operator  to  respond  as  required. 


10  OPTION  BASE  1  1  PROGRAM  FILE  “FLARE" 

lo  Din  a<2a, ;,R2<20>,*3 (io^,v<20> 

30  INTEGER  Arr«y  t T 0 : 23835 >, Array  o<0:23948> 

40  PRINT  PAGE 

50  BEEP 

80  PRINT  "THIS  PROGRAM  WILL  CALCULATE  AND  PLOT  ESTIMATED" 

70  PRINT  “  STATIC  AERODYNAMIC  COEFFICIENTS  OF  LONG  ROD" 

30  PRINT  “FLARED  PROJECTILES  IN  THE  MACH  RANGE  2-5" 

90  PRINT 


100  Pf*oJ_d*»i  <jns  PRINT 


1 10 
150 
130 
170 
130 
181 
190 
200 
210 
211 
220 
230 

240 

241 
250 
280 

270 

271 
280 
290 
300 
310 
320 
330 
340 
350 
430 
470 
480 
4?  I 

492 

493 
500 
510 

520 

521 
530 
540 

550 

551 

580 
570 
588 

581 
590 
308 

310 

311 
328 
338 


Coaput«d*0 

PRINT  "INPUT  THE  FOLLOWING  VALUES  IN  CALIBERS: “ 

PRINT  "  CONICAL  NOSE  LENGTH:  " 

BEEP 
INPUT  L 

PRINT  "L  *: “CONICAL  NOSE  LENGTH‘:L 
PRINT  "  CYLINDRICAL  BODY  LENGTH:  * 

BEEP 
INPUT  LI 

PRINT  "LI  "CYLINDRICAL  BODY  LENGTH'JLl 
PRINT  "  GROOVE  LOC  AT  I  ON ''NOSE:  " 

BEEP 
INPUT  PI 

PRINT  "PI  "; "GROOVE  L0CATI0n^N0SE“;P1 
PRINT  "  GROOVE  LENGTH;  * 

BEEP 

INPUT  P2 

PRINT  *P2  “GROOVE  LENGTH" ;P2 
PRINT  “  MAX  FLARE  DIAMETER" 

BEEP 

INPUT  H3 

PRINT  "H3  ":"MAX  FLARE  DIAMETER";H3 
PRINT  •  FLARE  ANGLE" 

BEEP 

INPUT  H4 

PRINT  "H4  "FLARE  ANGLE" :H4 

PRINT  “  CENTER  OF  GRAVITY  OF  PROJECTILE:  * 

BEEP 

INPUT  C2 

PRINT  "C2  "’.“CENTER  OF  GRAVITY  OF  PRO  JEC  T I LE  "  C  2 

PRINT 

PRINT 

PRINT  "INPUT  1.8  CALIBER  REFERENCE  DIAMETER  \ INCHES':  " 

BEEP 
INPUT  D 

PRINT  "D  "INPUT  l.O  CALIBER  REFERENCE  diameter  '.inches  '• ;  t 
PRINT  "INPUT  1.0  CAUSER  REFERENCE  DIAMETER  ■ MILLIMETERS: •  - 

BEEP 

INPUT  01 

PRINT  "  D 1  “{“INPUT  NORMALIZED  PROJECTILE  hEIC.HT  ■  M I LU  ME  TER  i  ' 
PRINT  "INPUT  NORMALIZED  PROJECTILE  WEIGHT  CUBIC  CALIBERS  •:  ' 

BEEP 

INPUT  U 

PRINT  “VI  ■;  "INPUT  NORMALIZED  PROJECTILE  WEIGHT  l.Sil  ‘HUBER 
PRINT  "INPUT  AXIAL  MOMENT  OF  INERTIA  . CAl I BEPS--T  ■ :  ' 

BEEP 
INPUT  I 

PPINT  "I  "{"INPUT  AXIAL  MOMENT  YF  I N£- "IP  C Al I ; E 1  I  - - T  111 
PRINT  "INPUT  TPHNSVEPSE  MOMENT  OF  INERTIA  C  Hl  I  BE  -  C  — 5  :  ' 

BEEP 


29 


^Vp%>Vc’S>'<v^'Sv.-’j"V  v-1' •/ 


64$  INPUT  II 

641  PRINT  “II  *' ;  “  INPUT  TRANSVERSE  MOMENT  OF  INEPTIA  ■  CAL  I  EEPS*»?  •  :  :  I  1 

650  PRINT 

660  PRINT  “INPUT  MACH  NUMBER  AT  MUZSLE:  " 

670  BEEP 
680  INPUT  M0 

631  PRINT  "M0  “INPUT  MACH  NUMBER  AT  MUSCLE: “ ; N0 
630  PRINT  “INPUT  MAXIMUM  RANGE  <  METERS > 1  “ 

700  BEEP 
710  INPUT  S 
720  S9»S 
730  PRINT  PAGE 
740  DEG 
750  B»20 

760  B 1 *25 
780  B3»L*L 1 
790  H6»H3-'2 

300  H7«1-'.5*(H3-1  > 

301  H8»B3+H7 

340  Image  l:  IMAGE  20.20 

850  2»0~ 

940  PLOTTER  IS  13, “GRAPHICS" 

950  lplot«0 

960  Graphics:  GRAPHICS 

970  EXIT  ALPHA 

990  SCALE  -6,40,5,23 

990  X  scat*  Tact  or»42  '23 

1000  MOVE  -1?20 

1010  ORAN  20,20 

1020  MOVE  0, B 

1030  ORAN  L.B+.S 

1040  DRAW  t.,  B-.  5 

1050  DRAU  0 , B 

1060  MOVE  L,B*.5 

1070  ORAN  B3.B+.5 

1030  ORAN  B3.B-.5 

1090  ORAN  L.B-.5 

1100  MOVE  B3,B*.5 

1110  ORAN  H7*B3.BrH6 

1120  ORAN  H7*B3,B-H6 

1130  ORAN  83.B-.5 

1540  MOVE  0,8-. 2 

1550  DRAM  0 , B-5 

1560  MOVE  L.B-.7 

1570  DRAW  L.B-2 

1530  MOVE  B3.B-.7 

1610  ORAN  83,8-4 

1640  MOVE  -. 25, B- 1 . 5 

1650  ORAN  L4.2S.B-i. 5 

1660  MOVE  -4, 8-1. 5 

1670  LORG  2 

1690  LABEL  L 

1690  MOVE  -.25,3-3 

1700  DRAW  834.25, B-3 

1710  IPLOT  -83-4. 3, 0,-2 

1720  LABEL  83 

1730  MOVE  -.25, B-3 

1740  OPAW  83*. 25. B-3 

1750  MOVE  -.25,8-4.5 

1763  DRAW  H9*. 25,5-4. 5 

1770  IPLOT  -H8-4. 3, O, -2 

1730  LABEL  H3 

1731  MOVE  H8. B-H6-.2 

1732  DRAW  H8.B-5 

1733  MOVE  H9.B-4.5 

1734  DRAW  n 3 , B - 4 . 5 


MOVE  0.B  +  .2 

DRAW  0,B  +  5 

MOVE  C2.B+3 

DRAW  C2.B+3 

MOVE  -.23,8+4.5 

DRAW  C2+.23.B+4.5 

IPLOT  -02-4.3,8, -2 

LABEL  02 

MOVE  -. 23, B+4. 3 

DRAW  02, B+4. 3 

MOVE  -. 23, 2  +  B 

DRAW  P1+P2, 2+B 

IPLOT  -P 1 -P2-4 . 3 , 3 , -2 

LABEL  PI 

MOVE  P1+P2, 2+B 

DRAW  P1+P2+.2.2+B 

MOVE  P 1 +P2+ . 3 , 2+B 

LABEL  P2 

LORG  1 

REM 

REM 

MOVE  C2-. 1, B-. 1 
DRAW  02+ . 1 , B+ .  12 
MOVE  C2+. 1.B-.12 
DRAW  02- . 1 , B+. 12 
MOVE  PI.B+.7 
DRAW  Pl.B+2.3 
MOVE  P1+P2.B+.7 
DRAW  P 1 +P2 , B+2. 3 
MOVE  02, B 
LORG  3 
LABEL  “0" 

MOVE  PI, .5+B 
DRAW  PI+P2, -. 3+B 
DRAW  P1+P2, .3+B 
DRAW  PI, -.3+B 
DRAW  PI, .3+B 
LORG  1 

MOVE  H8+. 23, B+. 3 

REM 

REM 

LABEL  H4 ; " DEG  FLARE" 

LORG  1 

MOVE  1.3.B-3 
LABEL  “WT  ■“ 

MOVE  4.3, B-9 
LABEL  WJ“  WCAL3” 

MOVE  1.3,8-13 
LABEL  "IX 
MOVE  4.3.B-10 
label  i;*  ical3 " 

MOVE  1.3.B-11 
LABEL  "IV  »“ 

MOVE  4.3,  B- 1 1 
LABEL  Ill"  ICAL3" 

MOVE  1.3,8-12 
LABEL  "  D I A » “ 

MOVE  4. 3,  B-12 

IF  D«0  THEN  LABEL  DU"  MM'  CAL" 

IF  D<  0  THEN  LABEL  Dl "  IN'CAL" 

MOVE  3.23,26 

LABEL  "LONG  ROD  FLARED  PROJECTILE  DESIGN" 
G3T0RE  Array _t.  ■'  +  ) 

GSTORE  Array  b' +  ',0,227 

EXIT  GRAPHICS 

ALPHA 


v‘v 


3330  IF  lplot*0  THEN  A»k_it 
3340  PLOTTER  7,5  13  OFF 
3330  GCLEAR 
3360  GOTO  Print_mput 
3370  A*k_it:  PRINT  PAGE 
3380  I p 1 ot  *  1 

3390  PRINT  "DO  YOU  WANT  A  HARD  COPY  OF  THE  PLOT'’  <Y  or  N>" 

3400  BEEP 

3410  INPUT  Reply* 

3420  IF  Rep  I y*< >“ Y“  THEN  Keep  picture 

3430  PRINT  “DO  YOU  WANT  THE  COPY  ON  THE  PRINTER  <0.',  OR  ON  THE  PEN  PLOTTER 

U 

3440  BEEP 

3450  INPUT  Plot_device 

3460  IF  Plot  deuice«0  THEN  Dump  gnpNici 

3470  PRINT  “SET  PEN  PLOTTER  AND_ PRESS  CONTINUE  WHEN  READY" 

3480  BEEP 
3490  PAUSE 

3300  PLOTTER  IS  ?,5,"?872A“ 

3310  GLOAt  Arr  »y_t  < *) 

3320  GLOAD  Arr4y~b«.*  J ,  0, 227 
3330  PLOTTER  13  Is  OFF 
3340  GOTO  Graphic* 

3330  Keep  picture:  PRINT  “PRESS  CONTINUE  TO  RETURN  PLOT.  THE  PLOT" 

3360  PRINT  "WILL  REMAIN  ON  THE  SCREEN  UNTIL  CONTINUE  IS  PRESSED  AGAIN." 

3370  BEEP 

3580  PAUSE 

3390  EXIT  ALPHA 

3600  GRAPHICS 

3610  PAUSE 

3620  GCLEAR 

3630  EXIT  GRAPHICS 

3640  ALPHA 

3630  GOTO  Print_mput 

3660  Dunp  graphic*:  EXIT  ALPHA 

3670  GRAPHICS 

3680  DUMP  GRAPHICS 

3690  GCLEAR 

3700  EXIT  GRAPHICS 

3710  Print  input:  PRINT  PAGE 

3720  PRINT  “DO  YOU  WANT  THE  INITIAL  DATA  PPINTED  OUT  “  <  Y  or  N>’ 

3730  BEEP 

3740  INPUT  Reply* 

3730  IF  Rep  I  >'*s  .  “  Y"  THEN  Menu 
3760  PRINTER  IS  0 
3770  PRINT 

3780  PRINT  "ALL  VALUES  ARE  IN  CALIBERS  UNLESS  OTHERWISE  NOTED" 

3790  PRINT 

3800  PRINT  “CONICAL  NOSE  LENGTH:  “;L 
3810  PRINT  “CYLINDRICAL  BODY  LENGTH:  “;L1 
3320  PRINT  “GROOVE  LOCATION^NOSE:  “;P1 
3830  PRINT  “GROOVE  LENGTH:  “ ; P2 
3840  PRINT  “MAX  FLAPE  D I AMETER: “ : H? 

3850  PRINT  "FLARE  ANGlE:  “;H4 

3900  PRINT  USING  If»«ge_l;F 

3930  PRINT  “CENTER  OF  GRAVITY:  ”;C2 

3940  PPINT  “MACH  NUMBER  AT  MU22LE:  " ; M0 

3930  PRINT  "MAXIMUM  RANGE  tMETERS):  ”;S 

3960  PRINT  “PROJECTILE  WEIGHT  c CALS':  “;W 

3970  PRINT  “AXIAL  MOMENT  OF  INERTIA  ■ CAL3  ■  :  “:I 

3980  PRINT  “TRANSVERSE  MOMENT  OF  INERTIA  ‘CALS’:  “;Il 

3990  IF  D-:  .‘0  THEN  °RINT  « i .  o  CALIBER  “EFEPENCE  DIAMETER  ‘In..;  •• :  D 

4000  IF  0 1  ,0  THEN  PRINT  "1.0  CALIBER  REFERENCE  DIAMETER  ■ KM  ■ :  " : D l 

4010  PRINT  L I N  <  3  > 

4020  PRINTER  13  16 
4030  Menu:  PPINT  PAGE 


4040  M tnu  list:  PRINT  “OPTIONS: “ 

4030  PRINT  “  0  -  END  PROGRAM* 

4060  PRINT  *  1  -  NOMENCLATURE* 

4070  PRINT  *  2  -  TABLE  OF  AEROBYNAM I C  COEFFICIENTS" 

4060  PRINT  “  3  -  CNAT  vs.  MACH  NUMBER  PLOT“ 

4090  PRINT  ■  4  -  CUT  vs.  MACH  NUMBER  PLOT" 

4100  PRINT  *  5  -  VELOCITY  us.  RANGE  PLOT* 

4110  PRINT  *  6  -  CMAT ,  CG  vs.  MACH  NUMBER  PLOT- 

4120  PRINT  *  7  -  J  ZETA  us.  MACH  NUMBER  PLOT- 

4130  PRINT  “  8  -  LONG  ROB  FLAREB  PROJECTILE  BESICN" 

4140  PRINT  "ENTER  SELECTION:* 

4130  BEEP 

4160  INPUT  Option 

4170  IF  Option»0  THEN  End 

4130  IF  Co«put«d<  >0  THEN  Uh*r**_to 

4190  Gl-U 

4200  G2-U 

4210  G6-L+L1 

4220  J2«L^S0R<L«L*.23) 

4230  J3«.3»<H3-1>'-TAN<H4> 

4240  M"2 

4230  FOR  Y-l  TO  7 
4260  G3-SQR<M*H-t> 

4270  C4-G3/-L 

4200  G3-L1/G3 

4290  Q1<Y)«M  _  .. 

4300  G?*<  1 . 9*1 . 3*G4* . 0149»G3>*( l'C3-.7>»< . 6 103*. 036*Ue 

4330  C8“TAN<F)/G3 

4340  G9“<  .  09*.  73*G4*.  3*G3>*<  l'G3'.S3>*<G3'M:'~2.3 

4330  J1-G7*G9 

4360  0<1,Y>«G? 

4370  Q(2,Y)"C9 

4300  Q<3,Y)»J1 

4390  J4*< 67/H3+C. 6?i»C0S<H4>*M^G3 

4400  J3*L*L1*.  1*<H3-1  )/"TAH<H4> 

4410  J6«J4*J3 

4491  J" J4*C7 
4300  Q<4,Y)«J4 

4310  K1«J 
4320  K2«I1*J6 

4330  K3»K2/'J 

4340  F4-C2-K3 

4330  K3*«4*J 

4360  K6«K3 

4370  Q<3,Y>«J 

4380  0<«,Y>"K2 

4390  0<7,y>»K4 

4600  Q<8,Y)*K3 

4610  J-6+J/10 

4620  K3»6*K3*6/220 

4630  K7»M*2-3.3 

4640  K8-J 
4630  K9«K? 

4660  Tl«K3 

4670  T2»2*L-' J2*4*L1  *4*  J  J-'COS  i  H4  >  ■♦<.  1  *-  J3*TANi  H4>  > 

4700  T3«H«.23»L~1.73 

4710  T3-.7/T3 

4711  T6»EXP< • 6*<  H4' 37. 3  >  > 

4720  T6«<.263-.043»M;«H3'2*T6 

4730  T7»T2*< . 084974- , 000721 #M> 

4731  T9«3*H4/37.3/M*<.7?-.6'H3> 

4780  F 1  ■»  00023*EXP  < 3 . 9*LQG<  M  > >*P2 
4790  F2*Fl*<T3*T7*T9> 

4300  F4«T3*T6*T7*T9*FI 

4318  P4«F4 

4820  REM 


4839 

4949 

4839 

4889 

4879 

4898 

4939 

4949 

4998 
4989 
4979 
4989 

4999 
3909 
3819 
3829 

3938 
3949 

3939 
3989 
3979 
3889 
3899 
3199 
3119 
3129 
3139 
3148  S 
3139 

3188 
3179 

3189 
3198 
3299 
3218 
3229 
3239 
3249 
3239 

3288 
3278 

3289 

3299 
3389 
3319 
5320 

3338 
3349 

3339 
3388 
3378 
3338 
3399 
3499 
3419 
3429 
3439 
3449 
3439 
t.S_p 
3480 
3479 
3*80 
3499 

3300 
3319 


REM 

Q<9, Y)*T3 

Q<10,Y)*T8 

0<11,Y>»T7 

0<12,Y>«F2 

9<14,Y)*T9 

9 '.18,  Y>*F4 

IF  N*3  THEN  F8-F4 

F4»F4»3*8 

IF  H*5  THEN  F7«<F4-8>'3 

F8»H*2*4 

F9-K1-F3 

8<19, Y)»F9 

N1«C2^C1»<F9-'KS1 

0<29,Y>*N1 

««H*.3 

Nl*8*Nl/2 

N2*Nl 

02-F8 

03-F4 

NEXT  Y 

N3-398 

Y-l 

N3»<F7-F6>'2 

N4»FS-3*N3 

N8-PI*.  87S//02.  4/3*N4/W 
N7»N3*N4/H0 

Skip  81  IF  B>9  THEN  N8»EXP< S*N8*39. 37/3) 

IF  31 >9  THEN  N8*£XP<S*H8*1808/31 > 

N9»N4/<N7*N8tN3> 

N9«<N9-N9) 434 1330 
V<Y)“N9*341 . 38 
N9-N0/S 
R2<Y)*$ 

R3<Y)-N8 
04-N9 
03-S 
S-S-N3 
Y-Y4  1 

IF  S<  >8  THEN  3kip_8 

N3-<F7-F6)/2 

Yl-Y 

R2<Y)-0 

R3<Y)*R3< Y-l ) 

V<Y)-H0*341 . 38 
N4*F8-3*N3 

N6-PI*. 073/62.  4/34N4  -H 
N7-N3*N4/M9 

IF  3)9  THEN  NS-EXP'OS'HSO-?.  37  D'  .*♦♦..*-♦4.-*-  1009" 

IF  31)0  THEN  N3-EXP< O3*N8«10O0' 31 ) 

N9*N4/ <N7*N9-N3) 

N8-<li0-N9)*341  380 
N0-N0/O3 

05« 1 1 /K6*8*P  t / .9012 
08»SQR <  MS  <  08 ) ) 

IF  B>0  THEN  01408*3/39.37 
IF  31 >0  THEN  01*08*31/1990 

Coaput  *d*  1  ( 

Hh#r*  toS  ON  Option  GOTO  Noa*nc  1  at  ur  » ,  Co«t'_ »  »b  '■  4 .  :  _P 1  o»  •  »_P 1  ,J'  •  •  »  _R  0 

plot.  Pro j 

S  plot:  'PRINT  PRGE 

PRINT  "30  YOU  Hi-iNT  Th£  PLOT  ON  THE  'SCREEN  0',  OR  ON  TrtE" 

PRINT  "  PEN  PLOTTER  <  1  > 

BEEP 

INPUT  Plo* 

IF  /Plot  9)  PHD  •  PI  o»  _1to  i  ■:  t  t1  TmEn  P '  o*  _  1* 1  ’  •:  *  »0 


3320  IF  Plot  d«u tc*»0  THEN  Skip  r*s«t 

3330  PPINT  "RESET  PLOTTER  PEN  LIMITS  FIND  PRESS  CONTINUE  WHEN  READY" 

3540  BEEP 
3330  PAUSE 

3360  PLOTTER  IS  7,3,‘9872A" 

3370  GOTO  Do  plot 

3380  Skip  r*s*  t:  PLOTTER  IS  13,  "GRAPHICS' 

3390  Do  plot:  GRAPHICS 
3800  EXIT  ALPHA 

5610  ON  Option  GOTO  Noa*nc 1 stur*, Co«f_t sbl *, Cnst , Cdt , 7# l_rsng» , Cast , J_:*t », Pro J 
_d«s1 gn 

3620  E  plot:  IF  Plot  d*uic*«0  THEN  End  scr**n 

5630  GCLEAR 

3640  PEN  0 

3630  EXIT  GRAPHICS 

3660  PLOTTER  7,5  IS  OFF 

3670  ALPHA 

3680  GOTO  N«nu 

3690  End  scr*«n:  EXIT  GRAPHICS 
3700  ALPHA 

5710  PRINT  "DO  YOU  WANT  A  COPY  OF  THE  PLOT  ON  THE  PRINTER ?  <Y  or  N)" 

3720  BEEP 
3730  INPUT  Copy* 

3740  IF  Copy*-"Y"  THEN  DUMP  GRAPHICS 

3730  PRINT  "PRESS  CONTINUE  TO  RETURN  PLOT.  THE  PLOT  WILL  REMAIN" 

3760  PRINT  "  ON  THE  SCREEN  UNTIL  CONTINUE  IS  PRESSED  AGAIN." 

3770  BEEP 

3780  PAUSE 

3790  EXIT  ALPHA 

3800  GRAPHICS 

3810  PAUSE 

3820  GCLEAR 

3830  EXIT  GRAPHICS 

3840  PLOTTER  13  IS  OFF 

3830  ALPHA 

3860  GOTO  Htnu 

3870  End:  PRINT  PAGE 

3880  END 


3890 

No«*ne 1 stur*: 

PRINT  PAGI 

3900 

PRINT 

"DO  YOU 

WANT  THE  1 

3910 

BEEP 

3920 

INPUT 

Printer 

3938 

IF  Print*r«0 

THEN  Do_ 

J  t 

3940 

PRINTER  IS  0 

3930 

PRINT 

L !N<  2 ) 

3960 

Do  it: 

PRINT  “ 

3970 

PRINT 

3980 

PRINT 

"CNAB 

SI  op* 

of 

dy* 

3990 

PRINT 

"XCPB 

Pr*s*ur* 

6000 

PRINT 

"CMAB 

S  1  op* 

of 

body 

i  ■ 

6010 

PRINT 

"CNAF 

Slop* 

OF 

sr  •  " 

6020 

PRINT 

"CNAT 

S!  op* 

Of 

6030 

PRINT 

"CHAT 

SI  op* 

of 

6040 

PRINT 

"CMAT.CG  Slop* 

of 

OR  ON  THE  PRINTER  '1 


NOMENCLATURE" 

Slop*  of  th*  Nor**!  Pore*  Co*ffici«n'  for  th«  proi*ctil*  bo 


:  1 1  1  *  b 
i  c i *nt 


odv" 
for  th* 


proj«o ’  i  I  « 

Slop*  OF  THE  Norast  Fore*  Co*ffi':i«nt  for  * h*  pro;*cti1«  ft 


6830 
6068 
6070 
6080 
6090 
6110 
6130 
6 1 60 
6170 


PRINT 
PRINT 
PRINT  “CDS 
PRINT  "CDV 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 


"CDUN 


"CDGR 
“CDWF 
"CDT 
"CLA 
"j  :eta 


c«nt*r  of  grsv  tty" 

Co*ffici*nt  of  usv*  drsg  for  th*  pr 
Sss*  drsg  eo*ffici*nt  for  th*  pro;* 
Viscous  dr sg  co*ffici«nt  for  th*  pr 
Profit*  drag  of  grcou*d  s*ct ion  of 
Wsv*  drsg  co*ff>ci*nt  for  U\«  proj* 
Tot  si  drsg  co*f  f  1  C  i  *nt  * 

Slop*  of  th*  lift  co*ffici*nt" 
A*rod'>nsf»  i  c  jump  fictor” 


»♦  f  i  ■:  i  • 
Co*t  f i 
Co*f f i 

o  ;*c »  1 
c  1 1  l  « “ 
o  ; « •: »  i 
bod'/ " 
c  ’■  1 1  * 


nt  “ 

c i *nt  “ 

ci*m  so 

i«  nos*" 
i*  “ 

f  1  sr*  " 


ou»  th* 


35 


»v» 

if* 

►fc: 

» * V 


i 

HP 

» *  j 


(>180  IF  Pr(nter»0  THEN  Beep 
6190  PRINT  LIN(53 
6288  PRINTER  IS  16 
6218  Beep:  BEEP 

6228  PRINT  LIN(13; “PRESS  CONTINUE  WHEN  READY  TO  RETURN  TO  MENU " 

6238  PAUSE 

6248  GOTO  H*nu 

6238  Coef  table:  PRINT  PAGE 

6268  PRINT  “DO  YOU  WANT  THE  TABLE  ON  THE  SCREEN  (03,  OR  ON  THE  PRINTER  ."l)?" 

6278  BEEP 

6288  INPUT  Printer 

6298  PRINT  PAGE 

6388  IF  Printer<>8  THEN  PRINTER  IS  8 

6318  PRINT  “  STATIC  AERODYNAMIC  COEFFICIENTS  FOR” 

6328  PRINT  <  LONG  ROD  FLAPED  PROJECTILES” 

6338  PRINT 

6348  PRINT  USING  laage  2; “MACH  NUMBER” 

6338  laage  2:  IMAGE  38x7llA 
6368  PRINT 

6378  PRINT  USING  laage  3; ”2. 8” , '2. 3“ , “3. 0“ , “3. 5“ . "4. 8” , “4. 3“ , “3. 0“ 

6398  laage  3:  IMAGE  12x73A,  6X,  3A,  7X,  3A,  7X,  4<  3A,  6X3 

6398  Iaage_4:  IMAGE  7A,  2(  IX,  4D.  3B3 , 2(2X,  4D.  3D3 , 3<  IX,  4D.  3D.. 

'  ”, 0(1,13,0(1, 23, Q(l, 33, 0(1,43.0(  1, 33, QM, 63, 0( 

•,Q(2,13,Q(2,23.Q'2,3),Q(2,4>,Q(2.3>,0(2,63.a< 

•,a(3,13,Q(3,23,Q(3,33,Q(3,4>,QC3,53,Q(3,«3,0' 

*,0(4, 13,G(4,23,Q(4,33,Q(4,43,Q(4,5>,Q(4,6 '.0' 

”,Q<3,  13,a(3,23,a(3.33,Q(5,43,Q<.3,3>,Q(5,6'<.Q( 

”,0<6,  13, 0(6, 23, 0(6,  3 ',0(6,  4>,  0(6,5  >,Qf  6.63, 0( 


6488 

PRINT 

USING 

I»»ge_4j “CNAB 

1.73 

6410 

PRINT 

USING 

I»*ge_4; ”XCPB 

2.73 

6420 

PRINT 

USING 

Iaage_4; "CNAB 

3,73 

6438 

PRINT 

USING 

Iaage_4; *CNAF 

4,73 

6440 

PRINT 

USING 

Ia*ge_4J "CNAT 

3,73 

6430 

PRINT 

USING 

Iaage_4; ”CMAT 

6,7) 

6468 

PRINT 

USING 

Iaage_4; *CG-CP 

7,73 

6470 

PRINT 

USING 

Iaage_4; *CMAT,CG 

8,73 

6480 

PRINT 

USING 

I»*ge_4; "CDWN 

9,73 

6498 

PRINT 

USING 

I«4ge_4; ”CDB 

,63 , 

10(10,73 

6300 

PRINT 

USING 

I«*ge_4; “CDV 

.63, 

0<U,73 

6310 

PRINT 

USING 

Iaage_4; "CDGR 

.63, 

0(12,73 

6330 

PRINT 

USING 

I»age_4; "CDWF 

•  6  J , 

0(14,73 

6370 

PRINT 

USING 

Ia*ge_4; “CDT 

.63, 

0(19,73 

6389 

PRINT 

USING 

I»*ge_4: ”CLA 

.63, 

0(19,73 

6398 

PRINT 

USING 

Iaage_4;”j  2ETA 

,6), 

0(28,73 

',0(9,  13, 0<9,2  3,0(9,  3  3,  Q(»,  4  3, 0(9,3  3,0  (9,63,  Q( 

”,  0(  18,  1  >,0<  18.2',3'  10,  3>  .  0<  13, 4  i ,  0>  10, 5'-  ,0'  10 

”,0<  It ,  13  ,0<l  l,23,0(  11 . J3,Q(  1 1 . 4  ■  .0-.  1 1 . 3  >,  0-  1 1 

',0(12, 1 3.0< 12,23,0< 12, 33 ,0( 12,4 ' , 0.  12,3  '  .0'  12 

",Q(14, 1 3,0( 14,23,0( 14,33 ,Q( 14, 4 ',0< 14, S> ,0: 14 

”,Q( 18, 1 >,0( 18,23,0*  18,33,0' 13, 4  ,0'  16, 3  ..0.  13 

1 ,  Q>  19,  l  >  ,0(  19. 2  >,00  9, 3>  .O'.  19, 4  •  .0(1  9.3  • .  O'  19 

•J  2ETA  ",0(28,  l;, 0*20, 2 ',0(20.3 '.0*20,4  ,0<  20,5.  .0' 20 

6608  IF  Printer»0  THEN  Skip  7 
6618  PRINT  LIN<33 
6628  PRINTER  IS  16 
6638  GOTO  Menu 
6648  Skip  7:  BEEP 

6638  PRINT  LIN<13;  "PRESS  CONTINUE  WHEN  READY  TO  RETURN  TO  MENU" 

6668  PAUSE 
6678  GOTO  Menu 
6690  Cn«t :  SCALE  1 . 4 , 3. 2, -6, 3 1 
6690  CLIP  2,3,9,38 
6788  LAXES  .23,2,2.8, -4, 3 
6718  UNCLIP 
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6720  LOIR  0 

6730  LORG  5 

6740  HOVE  3. 3, -3 

6730  LABEL  "MACH  HUMBER “ 

6760  LOIR  90 
6770  HOVE  1.7,13 
6780  LABEL  "CHAT" 

6790  MOVE  fll< 1),0<3, t> 

6800  FOR  1-2  TO  7 
6810  9RAW  Q1 ( I ) , Q(3,  I) 

6820  NEXT  I 

6830  GOTO  E  plot 

6840  Cdt:  SCALE  1 . 4, 3. 2, 2, 1 . 03 

6830  CLIP  2,3,0, 1 

6860  LAXE3  . 23, . 1 , 2, 0, -4, 2 

6870  UNCLIP 

6880  LOIR  0 

6890  LORG  3 

6900  MOVE  3.3, -.1 

6918  LABEL  “MACH  HUMBER" 

6920  LOIR  90 
6930  MOVE  1.7,. 3 
6940  LABEL  "CDT" 

6930  MOVE  Ql<l),Qtl8,  1> 

6960  FOR  1-2  TO  7 
6970  DRAU  01 < I > , 0< 18, I > 

6980  HEXT  1 
6990  GOTO  E_j»1  Ot 

7000  V*1  rang*:  SCALE  -S9x3000, 3 l *S9'30O00, - 1600, 2400 

7oie  clip  e,s9/ieee, 0,2000 

7020  LAXES  .3,230,0,0,-1,  1 

7030  UNCLIP 

7040  L9IR  0 

7030  LORG  3 

7060  MOVE  S9'2000,-23e 

7070  LABEL  "RANGE  < KM) “ 

7080  LSIR  90 

7090  MOVE  -S9/-8000,  1000 

7100  LABEL  "VEL  <M'S>" 

7110  HOVE  R2< 1 >'1000, V< 1 ) 

7120  FOR  1-2  TO  Y1 
7130  DRAU  R2U  >'1080,  V<  I  > 

7140  NEXT  I 
7130  LOIR  0 
7160  LORG  1 
7170  MOVE  0,2130 

7ia0  Ia*g*  3!  IMAGE  “MUZZLE  MACH  No . ” , 2D . 30 , " ,  INITIAL  VAN  CYCLE" 
7190  LABEL  USING  la*?*  3*,H0,01 
7200  MOVE  0,-300 

7210  Ia«g*  6:  IMAGE  "RANGE", 3D,"  M,  RETARDATION" , 4B. 3D, "  M5-KM" 

7220  FOR  I-Yl  TO  1  STEP  -l 

7230  LABEL  USING  !**g*  6; R2 <  l  > , R3<  I ) 

7240  IPLOT  0,-30, -2 

7230  NEXT  I 

7260  GOTO  E  plot 

7270  Ca«tl  SCALE  1.4,3.2,10,-23 

7280  CLIP  2, 3, 6, -22 

7290  LAXES  .23, 2, 2, 0,-4, 2 

7300  UNCLIP 

7310  LOIR  0 

7320  LORG  3 

7330  MOVE  3.3,3 

7340  LABEL  "MACH  NUMBER" 

7330  LOIR  90 
7360  MOVE  1. 33,-0 
7370  LABEL  "CMAT.CG" 


4D. 3D, "M" 
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r-- 

7460 

MOVE  Q1<1>,Q<3, t> 

7470 

FOR  1 -2  TO  7 

7480 

ORAM  Q 1 < I > , Q<  8, I) 

|V 

7490 

NEXT  1 

[K 

7900 

GOTO  E  plot 

7910 

J  z»t SCFH-E  1.4, 9. 2, -2, 

■fl 

7920 

CLIP  2,9,0,20 

■ 

7930 

LAXES  .29, 1,2,0, -4, 2 

7940 

UNCLZP 

L> 

7990 

LOIR  0 

H> 

7900 

LORG  9 

a) 

7970 

MOVE  3.9, -1 

w 

7980 

LABEL  ’MACH  NUMBER’ 

Ul 

7990 

LOIR  90 

■ 

7000 

MOVE  1.0,10 

1 

7010 

LABEL  *J  2ETA* 

PP 

7020 

X  Inc *1.8 

By 

7030 

LOIR  8 

rj 

7040 

LORG  2 

fry 

7090 

FOR  1-2  TO  20  STEP  2 

lV 

7000 

MOVE  X  me,  I 

[N 

• 

7070 

LABEL  T-‘ 

fv 

7080 

IF  I >7  THEN  X  mc-1.79 

U 

7090 

NEXT  I 

9 

7700 

MOVE  Q1<1),ABS(Q<20, 1>> 

7710 

FOR  1-2  TO  7 

7720 

ORAM  Ql(D,ABS<Q<20t  I>> 

7730 

NEXT  I 

& 

7740 

GOTO  E _plot 
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